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We present measurements of the branching fractions for the decays B ± — > af (1260)7r° 
and B ± — > 05(1260)^* from a data sample of 232 x 10 6 BB pairs produced in 
e + e~ annihilation through the T(4S) resonance. We measure the branching fraction 
B(B ± -> of (1260)7r°) x B(af (1260) -> 7r"7r+7r ± )=(13.2 ± 2.7 ± 2.1) x 1(T 6 with a significance of 
4.2cr, and the branching fraction B(B ± -> a?(1260)7T ± ) x B(o?(1260) -> 7r"7r + 7r ) = (20.4 ± 4.7 ± 
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3.4) x 10 6 with a significance of 3.8a, where the first error quoted is statistical and the second is 
systematic. 



PACS numbers: 13.25.Hw, 12.39.St, 11.30.Er 

The rare decays of B mesons to two-body final states 
with an a!(1260) and a ,tt° , K^ 1 or K° are important 
processes for testing theoretical factorization model pre- 
dictions for branching fractions, branching fraction ratios 
and CP-violation parameters. The measurements can be 
combined with assumptions about SU(3) symmetries to 
form upper bounds on Aa =| a — a e ff | , where a is the 
weak interaction phase a = arg [—V td V t * b /V ud V* b ] of the 
Unitarity Triangle [l| and a e ff is the measured phase. 
The difference Aa is a measurement of the poorly known 
strength of the penguin amplitudes in the decay and can 
be used to improve our understanding of the CP- violating 
mechanism. 

The rare decays B ± -> 0^(1260)71° and 
B ± -> a?(1260)7r ± are expected to be dominated 
by b — > uud contributions. The branching frac- 
tion for B° — > af 7r T has been measured to be 
(33.2 ±3.8 ±3.0) x 10" 6 0] and this agrees well 
with the calculation of Bauer, Stech and Wirbel [3} 
within the framework of naive factorization and as- 
suming iKib/Kbl = 0.08. A more recent analysis 
using naive factorization and measured form factors 
predicts branching fractions in the range (5 — 11) x 10~ 6 
and (4 - 9) x 10~ 6 for B ± -► a^Tr and B ± -> a?7r ± , 
respectively jfj}. Previous measurements have placed 
90% confidence level upper limits of 1.7 x 10~ 3 and 
9 x 10~ 4 on the branching fractions for B ± — * 0^71° and 
P ± — > a°7r ± , respectively [Bj], and recently the BABAR 
collaboration reported the first measurements of the 
CP-violating asymmetries in the decay B a — > 7r T 0] . 

We present measurements of the branching fractions 
for the two charmless B meson decays P ± — ► aftr and 
B^ — > a^iT^ where the final state contains one neu- 
tral and three charged pions. The a\ — > 37r decay pro- 
ceeds mainly through the intermediate states (7r7r) p 7T and 
(7r7r) CT 7r Q. We do not distinguish between the dominant 
P-wave (7T7t)p and the S-wave (jnr) a in the channel 7r + 7r~. 
Possible background contributions from B — > a2(1320)7r 
are investigated. Charge conjugate modes are implied 
throughout this paper. 

The data were collected with the BABAR detector Q 
at the PEP-II asymmetric e + e~ collider. An integrated 
luminosity of 211 fb , corresponding to 232 million 
BB pairs, was recorded at the T(AS) resonance ("on- 
resonance") at a center-of-mass (CM) energy y/s = 
10.58 GeV. An additional 20 fb -1 were taken about 
40 MeV below this energy ( "off-resonance" ) for the study 
of continuum background in which a charm or lighter 
quark pair is produced. 

Charged particles are detected and their momenta 
measured by the combination of a silicon vertex tracker, 



consisting of five layers of double-sided silicon detectors, 
and a 40-layer central drift chamber, both operating in 
the 1.5-T magnetic field of a solenoid. The tracking 
system covers 92% of the solid angle in the CM frame. 
Charged-particle identification (PID) is provided by the 
average energy loss (dE/dx) in the tracking devices and 
by an internally reflecting ring-imaging Cherenkov de- 
tector. A K/ir separation of better than four standard 
deviations (a) is achieved for momenta below 3 GeV/c, 
decreasing to 2.5a at the highest momenta in the B decay 
final states. 

The off-resonance data together with the Monte Carlo 
(MC) simulations of the signal decay modes, continuum, 
BB backgrounds and detector response [3] are used to 
establish the event selection criteria and reconstruction 
efficiency. The MC signal events are simulated as B + 
decays to ai7r with a\ — > pn. The a\ and 02 line shapes 
are generated with EvtGen [ic|. where we use mass and 
width parameters from Refs. [2| and [7J. 

Two photons with a minimum energy of 30 MeV 
(100 MeV for B + — > a\ir + ) and an invariant mass of 
120 < m 77 < 150 MeV/c 2 are used to reconstruct the 
7r°. The intermediate dipion states (7r + 7r~) or (tt+tt ) 
are required to have an invariant mass of 0.46 < m^n < 
1.1 GeV/c 2 . We impose PID requirements to cleanly 
identify the charged pions and to suppress contamination 
from a\K. We require the invariant mass reconstructed 
for candidate — ► 7r~7r + 7r + and a\ — ► 7r~7r + 7r° decays 
to be 0.8 < m ai < 1.8 GeV/c 2 . 

A B meson candidate is characterized kine- 
matically by the energy-substituted mass toes = 
\J (a/2 + p • pb) 2 /Eq — p 2 B and energy difference 
AE = E* B — i/s/2, where the subscripts and B refer 
to the initial T(4S) and to the B candidate in the 
lab-frame, respectively, and the asterisk denotes the 
Y(AS) frame. The resolutions in toes an d in AE are 
about 3.0 MeV/c 2 and 20 MeV, respectively. Candidates 
are required to have 5.25 < toes < 5.29 GeV/c 2 and 
\AE\ < 0.2 GeV. To reduce fake B meson candidates we 
require a B vertex x 2 probability > 0.01. The absolute 
value of the cosine of the angle between the direction 
of the 7T meson from a\ —* pn with respect to the flight 
direction of the B in the a\ meson rest frame is required 
to be less than 0.85 to suppress misreconstructed candi- 
dates. The distribution of this variable is flat for signal 
and peaks near ±1 for misreconstructed candidates. 

To reject continuum background, we use the angle 6t 
between the thrust axis of the B candidate's decay prod- 
ucts and that of the rest of the tracks and neutral clusters 
in the event, calculated in the CM frame. The distribu- 
tion of cos#t is sharply peaked near ±1 for combinations 
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drawn from jetlike qq pairs and is nearly uniform for the 
isotropic B meson decays; we require | cos#t| < 0.65. 

The decay mode B — > a 2 7r can also give background 
contributions. It is suppressed by using the angular vari- 
able A, defined as the cosine of the angle between the 
normal to the plane of the 3tt resonance and the flight 
direction of the bachelor pion evaluated in the 37r reso- 
nance rest frame. Since the a± and a 2 have spins of 1 
and 2, respectively, the distributions of A for these two 
resonances differ. We require \A\ < 0.6, which reduces 
the a 2 background by more than a factor of two in both 
decay channels. 

After all the above selections, we have on average 1.20 
and 1.56 candidates per event in events where there is 
at least one candidate, for B + — > a^n° and B + — > a\ir + , 
respectively, and we select the B candidate with the (7777) 
mass nearest to the nominal p mass [7]. From the sim- 
ulation, we find that this algorithm selects the correct- 
combination candidate in B + — > a~l it and B + — ► a\-K + 
in 65% and 55% of events containing multiple candidates, 
respectively. 

We use an unbinned maximum-likelihood fit using five 
variables to extract the background and signal yields of 
B + — > a+7r° and B + — > a\i: + . We describe the B decay 
kinematics with the two variables AE and toes- We also 
include the invariant mass of the 37r system (m ai ), the 
variable A and a Fisher discriminant T . This discrimi- 
nant combines four variables: the angles with respect to 
the beam axis of the B momentum and B thrust axis in 
the CM frame, and the zeroth and second angular mo- 
ments of the energy flow around the B thrust axis Q ■ 

The extended likelihood function is 



N 



£= iV! exp 



n 



3 



flj"Pj(Xi) Otj) 



(1) 



where rij is the yield of events for hypothesis j (signal, a 2 , 
BB charmless, BB charm or continuum) and N is the 
number of events in the sample. The probabilities Vj are 
products of probability density functions (PDF) for each 
of the independent variables Xi — {iwes, AE, m ai , T ', A} 
evaluated for each event i. The otj are the parameters 
of the distributions in By minimizing the quantity 
— In C in two separate fits, we determine the yields for 
B+ -> a+7r° and B+ -> a?7r+. 

To take into account the relatively large number of mis- 
reconstructed signal events, the signal is separated into 
two components, representing the correctly reconstructed 
(true) and the self cross-feed (SCF) candidates, with pro- 
portions fixed in the fit for each mode. SCF occurs when 
a track from an cl^tt or a^Tr" 1 " is exchanged with a track 
from the rest of the event. The fraction of SCF, deter- 
mined from MC, is 35% and 44% for B+ -> tt° and 
B + — > a]*7r + , respectively. 

In addition to the a 2 , there are three main categories 
of backgrounds: BB charmless, BB charm and contin- 



uum. BB backgrounds are studied using MC simula- 
tions of B a B° and B + B~ decays, using a large sample 
equivalent to ~ 0.8 ab -1 . There are 17 BB charmless 
decays for B + — ► tt° and 20 for B + — > a\ir + that con- 
tribute as background. Those decays with similar distri- 
butions are grouped to form 13 and 10 hypotheses, re- 
spectively, and are included in the fit with a fixed yield as 
determined from MC. The total BB charmless yields are 
368 ±92 and 755 ±164 for B+ -> a^7r° and B+ -> a?7r+, 
respectively. These are dominated by B — > pp, B — > a±p 
and the other B — ► aiir mode under study. The BB 
charm backgrounds are included as a single hypothesis, 
with the normalization of the BB charm yield as a free 
parameter. Continuum events come from light quark pro- 
duction. We establish the functional forms and parame- 
ter values of the PDFs for BB charm and BB charmless 
backgrounds from MC simulations. For continuum, we 
use off-resonance data for the Fisher, on-resonance data 
with \AE\ > 0.1 GeV for toes, and on-resonance data 
with 5.25 < tties < 5.27 GeV/ c 2 for the other variables. 

We model the distributions using appropriate func- 
tions. The A distributions are modeled with polynomials. 
For the true signal component, the remaining distribu- 
tions are fitted using modified Gaussians ll| , and a rel- 
ativistic Breit-Wigner line-shape with a mass-dependent 
width [12], as necessary. The SCF component and the a 2 
have similar shapes to the true signal but have broader or 
more asymmetric distributions and shifted means. The 
BB backgrounds and continuum distributions are mod- 
eled with modified Gaussians, polynomials, nonparamet- 
ric functions 13j and, for toes , a phase-space- motivated 
empirical function [14|. The PDF variables are assumed 
to be independent except for B + — » a\'K + , where a two 
dimensional nonparametric PDF 13] in m ai and AE ac- 
counts for observed correlations in the MC for both true 
signal events and SCF. 

In the fit there are six free parameters: four yields 
(signal, continuum, a 2 and BB charm background), and 
two continuum background parameters (AE polynomial 
coefficient and toes shape coefficient £ [14| ). 

For B + — > a^TT , there are 24608 events in the data 
sample. We measure the raw signal yield to be 459 ± 78 
events with a reconstruction efficiency of 12.5 ± 0.1%, 
corrected for differences in tracking and neutral particle 
reconstruction between data and MC. The yield of the 
decay B + -> a^ir is 28 ± 65 events. For B+ -> 0?^+, 
there are 33375 events in the data sample and we measure 
the raw signal yield to be 382 ±79 events with a corrected 
reconstruction efficiency of 7.2 ± 0.1%. The yield of the 
decay B + — * aS L it + is 107 ± 65 events. 

We confirm our fitting procedure by generating and fit- 
ting MC samples containing signal and background pop- 
ulations using the yields as found from data. We iden- 



tify a signal yield bias for £T 



and B' 



of 16.8 ± 0.1% and 10.9 ± 0.1%, respectively. We fit 
for the branching fractions taking into account the fit- 
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ted signal yield, the yield bias, the corrected reconstruc- 
tion efficiency, daughter branching fractions, and the 
number of produced B mesons, assuming equal produc- 
tion rates of B°B° and B + B~ pairs. The statistical 
significance is taken as the square root of the differ- 
ence between the value of — 21n£ for zero signal and 
the value at its minimum. We measure the branching 
fraction B(B+ -► o+tt ) x B(af -► 7r-7r+7r+) = (13.2 ± 
2.7) x 10~ 6 with a statistical significance of 5.3cr and the 
branching fraction B(B + — ► a"-7r + ) x B{a\ — ► n~n + n ) 
= (20.4 ± 4.7) x 10~ 6 with a statistical significance of 
4.7cr, where the errors are statistical. 




1.4 1.6 1.8 

m ai (GeV/c 2 ) 



Fisher 



Projections of a) AE, b) tties, c) 



FIG. 1: (color online). 

m ai , and d) T for B + — » a^iv 1 '. Points represent on-resonance 
data, dashed lines the signal, dotted lines the continuum, 
dashed-dotted lines the BB charm background, the filled re- 
gion the ei2 background and solid lines the full fit function. 
These plots are made with a requirement on the signal likeli- 
hood to enhance the signal, and thus do not show all events 
in the data sample. 

Figs. [T]and [2] show the AE, ttt-es, m ai , and T projec- 
tions for B + — > a~l 7T° and B + — > o!{tt + made by selecting 
events with a signal likelihood (computed without the 
variable shown in the figure) exceeding a threshold that 
optimizes the expected sensitivity. 

The systematic errors are summarized in Table H] We 
determine the sensitivity to the parameters of the sig- 
nal and background PDF components by varying these 
within their uncertainties. The effect of varying the mass 
and width of the eii by the errors as reported in Rcf. [2] 
is included in the PDF parameters' variation systematic. 
The uncertainty in the fit bias correction is taken as half 
of the fit bias correction. The effect of possible inter- 
ference between ai and a\ is estimated by adding the 



5.27 5.28 5.29 
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FIG. 2: (color online). Projections of a) AE, b) tues, c) m ai , 
and d) T for B + — > a?7r + , using the same criteria and line 
styles as Fig. \T\ 



TABLE I: Summary of systematic errors for the a^Tv and 
branching fraction measurements. 



Systematic 


a+7r U 


aV+ 


PDF Parameter Variation 


8.6% 


8.8% 


Fit Bias 


8.4% 


5.5% 


£8i — a,2 Interference 


6.6% 


7.4% 


SCF Variation 


4.4% 


8.2% 


Tracking Efficiency 


3.9% 


3.9% 


7T° Efficiency 


3.0% 


3.0% 


Flight Direction Criteria 


2.0% 


2.0% 


P-wave and S-wave Reconstruction 


1.6% 




Charmless BB Background 


1.4% 


3.1% 


Number of BB Pairs 


1.1% 


1.1% 


cos 6t Selection Criteria 


1.1% 


1.8% 


Track Multiplicity 


1.0% 


1.0% 


pnn, 4tv Cross- Feed 


0.9% 


0.5% 


aiK Cross-Feed 




0.4% 


Total 


16% 


16% 



ci2 and ai amplitudes together with a varying phase dif- 
ference and using half the maximum change in yield as 
an uncertainty. The uncertainty in SCF is investigated 
by varying the SCF fraction. We also perform a sepa- 
rate fit treating the SCF as an independent background 
component. The fitted branching fraction is compati- 
ble with the nominal fit within the increased statistical 
uncertainty, but the statistical significance is reduced to 
3.5a and 3. Oct for B + — * a^n° and B + — * a\it + , respec- 
tively. A systematic uncertainty of 1.6% is estimated for 
the difference in reconstruction efficiency in the decay 
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modes through the dominant P-wave (nir) p and the S- 
wave (tttt),^ An error is assigned for the uncertainty in 
the fixed charmless BB background yields and possible 
interference effects by varying the individual components 
by the reported error on the branching fractions [7[ • The 
systematic errors for the flight direction criteria, number 
of BB pairs, cos#t selection criteria, track multiplicity, 
potential backgrounds from pirir and 47r, and a\K cross- 
feed are small. The total systematic error for both modes 
is 16%. The significance of the branching fractions, com- 
bining both statistical and systematic errors, is 4.2a for 
B+ -> a+7r° and 3.8cr for B+ -> a?7r + - 

In conclusion, we have measured the branching frac- 
tions B(B ± of (1260)tt°) x B(af (1260) -► tt-tt+tt ± ) 
= (13.2±2.7±2.1) x 10~ 6 and B(B ± -> a?(1260)7r ± ) x 
B(a?(1260) -> tt-tt+tt ) = (20.4 ± 4.7 ± 3.4) x 10~ 6 . 
Neglecting isoscalar contributions to the two-pion 
state, we assume B(a x (1260) — > 7r _ 7r + 7r ± ) is equal to 
B(af (1260) -> t^tt ) and B(af (1260) -> (3tt) ± ) is 
equal to 100% 0, resulting in 5(5="= -► of (1260)tt°) 
= (26.4 ± 5.4 ± 4.1) x 10~ 6 . We measure 
B{B ± -> a?(1260)7r ± ) = (20.4 ± 4.7 ± 3.4) x 10~ 6 , 
assuming B(o?(1260) -> tt-tt+tt ) is equal to 100%. The 
first errors quoted are statistical and the second are 
systematic. The signals are seen with significances of 
4.2c and 3.8a, respectively, and are in agreement with 
factorization model predictions 
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